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Abstract 

We explore the potential of several Neutrino Factory (NF) setups to constrain, discover and 
measure new physics effects due to Non-Standard Interactions (NSI) in propagation through Earth 
matter. We first study the impact of NSI in the measurement of 613: we find that these could be 
large due to strong correlations of ^13 with NSI parameters in the golden channel, and the inclusion 
of a detector at the magic baseline is crucial in order to reduce them as much as possible. We 
present, then, the sensitivity of the considered NF setups to the NSI parameters, paying special 
attention to correlations arising between them and the standard oscillation parameters, when all 
NSI parameters are introduced at once. Off-diagonal NSI parameters could be tested down to the 
level of 10"'^, whereas the diagonal combinations (eee — Crr) and (e^^ — e^r) can be tested down 
to 10^^ and 10^^, respectively. The possibilities of observing CP violation in this context are also 
explored, by presenting a first scan of the CP discovery potential of the NF setups to the phases 
(peii,(peT and 6. We study separately the case where CP violation comes only from non-standard 
sources, and the case where it is entangled with the standard source, 6. In case 6 turns out to be 
CP conserving, the interesting possibility of observing CP violation for reasonably small values of 
the NSI parameters emerges. 

PACS numbers: 14.60.Pq,14.60.Lm 
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I. INTRODUCTION 



From the results of neutrino oscillation experiments jlH3| we now know that neutrinos 
have masses and there is flavor mixing in the lepton sector [J] . However, the leptonic mixing 
is still not completely understood. In the three-family u-mass enriched Standard Model 
(i^SM), there are still three unknown oscillation parameters: 6'i3, the CP- violating phase 
5 and the mass hierarchy sgn{Am1j). Several reactor [5|-0] and accelerator jsi [9| neutrino 
experiments are currently running, or will start running in the near future, to search for 
positive signals of non-zero ^13. Strategy for exploration of the remaining two unknown 
parameters in the leptonic mixing, 6 and the mass hierarchy, heavily depends upon whether 
they succeed or fail to detect non-zero 6*13. 

If sin^26'i3 < 10~^, the ongoing and the near future experiments would face difficulties 
to observe clear evidence for 61^. In this case, a new generation of experiments will be 
needed, beyond doubt, to determine the remaining unknowns in the leptonic mixing. The 



candidates for such facilities are the High- Energy Neutrino Factory (HENF) [10|, the 
Beta-Bea m |l2i[l3t . and upgraded conventional beams {a.k.a. SuperBeams; see, for example, 
Refs. 0, [IITIiir The main motivation for a facility as ambitious and technologically 
demanding as the HENF is commonly considered to be the search for a very small ^13 
(sin^ 26*13 ^ 10""^) [l6[ and, consequently, the possibility to measure S and the mass hierarchy 
in this regime, something beyond the reach of both Beta-Beams and Super-Beams (l7| . 

On the contrary, if ^13 is relatively large so that the forthcoming experiments are able 
to see its effects, alternative strategies and scenarios for further exploration of the leptonic 
mixing could be possible. Such a possibihty of 6*13 being of the order of the Chooz limit 
18-21] has been recently suggested by some global analyses [22-28|| but it is still controver- 
sial. If confirmed, possible strategy and scenarios for further exploration of lepton mixing 
would have to be changed. A possibility is that facilities such as upgraded SuperBeams or 
Beta-Beams could be exploited to measure the standard oscillation parameters. Another 
possibility is that different designs for the Neutrino Factory (NF) can be exploited in or- 
der to achieve the same goal. Such a scheme is, for example, represented by the so-called 
Low-Energy Neutrino Factory j29|, 30] . However, the physics case for a HENF scheme can 
be well motivated even in the case of a relatively large 6*13, given its enormous accuracy. In 
particular, the HENF could be re-designed in order to look for both the i^SM parameters 
and possible New Physics (NP) beyond z/SM. These include, for example, Non-Standard 
neutrino Interactions (NSI) j31-35||, effects of non-unitarity of the leptonic mixing matrix 
3 [13, and new phenomena due to light sterile neutrinos [38[. Notice that, for a relatively 



small ^13 (sin^ 26*13 ^ 10~^), the NP parameters discussed above can spoil the sensitivity of 
the HENF to standard oscillation parameters, a possibility that must be studied. Recent 
works dedicated to the study of the potential of HENF in this context include, for example, 



Refs. [39|-|45[. For an extensive list of remaining references see, e.g., Ref. [46 . 

In the present work, we will attempt a complete study of the NSI effects in neutrino 
propagation. Since this kind of NSI can be regarded as a new "effective matter potential" , a 
HENF with very long baselines seems to be the optimal facility to study their effects. Due 
to an advanced minimum-searching algorithm in multi-parameter space, the MonteCUBES 
software [4^, we are able to include, for the first time, all NSI parameters in propagation (in- 
cluding their associated CP-violating phases) at the same time into our analysis. Therefore, 
we will pay special attention to correlations appearing when the whole set of NSI parameters 
is introduced at once in the simulations, a problem that has never been addressed in the 
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literature. First, we will study the potential effect of these NSI in the measurement of the 
standard parameters ^13 and 6. Then, we will examine the sensitivity to the NSI parameters, 
studying their correlation with the whole set of parameters, standard and NSI ones. Finally, 
we will explore the new avenues of CP violation coming from NSI and the relation among 
the different CP-phases involved. 

Intensive efforts have been devoted to find out which of the proposed future neutrino 
oscillation facilities is the best option. The International Scoping Study (ISS) ended with 
the publication of three reports on the status of HENF and their competitors in 2007 fv^ 



481 . 14 9| . Currently, the study is being continued by the International Design Study for a 
Neutrino Factory (IDS) [50], which has proposed a "baseline" HENF setup which optimizes 
the performance of the NF for measuring the uSM parameters: it consists of a neutrino 
beam obtained from the decay of 25 GeV muons aimed at two magnetized iron neutrino 
detectors (MIND) located at 4000 km and 7500 km. This setup makes use of two channels: 
the "golden channel" 16|, z/g — > u^, and the muon disappearance channel z/^ — )■ z/^, together 
with their CP conjugate ones. We will study the sensitivities than can be achieved at this 
particular setup in the context of NSI. 

As stressed before, however, new HENF designs can be explored to optimize searches 
for NP beyond the z/SM parameters. In particular, hi ghe r neutrino energies turn out to be 
better in order to study effects of NSI in propagation [40| , since these can be regarded as a 
kind of generalized matter effect. We will also study, therefore, a simple modification of the 
setup described above by increasing the parent muon energy to 50 GeV. In the context of 



NSI sensitivity studies, a quite similar setup was previously examined in Refs. |40l443 . 

In this paper, we also study a different setup with a composite detector consisting of 
a MIND and a Magnetized Emulsion Cloud Chamber (MECC) located at 4000 km from 
the source, with doubled statistics with respect to IDS-inspired proposals in which each of 
the two baselines receives half of the available neutrino flux. The MECC is an enlarged and 
more sophisticated version of the emulsion detector used in the ongoing OPERA experiment 
[5l| . which aims at detecting r particles with much higher efficiency. Using this setup, it is 
possible to measure the NSI parameters using two additional channels: Ur, the "silver 

channel" 52|, and z/^ — ?■ Ur, the "discovery channel" 45|. In order to overcome the strong 



suppression due to the small z/^-A^ cross-section, a parent muon energy of 50 GeV is adopted 
also in this case. Notice that for z/SM parameter searches, it was shown that this setup is 



outperformed by the baseline IDS setup defined above (see Refs. [17|, l53[). In Ref. |4l| this 
analysis was extended to the case of NSI in propagation when one NSI parameter is turned 
on at a time, with similar results (see, also, Ref. 54| for a summary on the performance of 
this setup in these cases). We want to check in this paper if these results hold in the case in 
which all NSI parameters are turned on simultaneously, i.e. in a case in which correlations 
between z/SM and NSI parameters are taken into account. 

The paper is organized as follows. In section [Tll the formalism for NSI is presented, and 
the main dependences of the probabilities on the parameters are introduced; in section IIIII 
we introduce the statistical approach we have used, and the details for the three NF-based 
setups we have studied; section [IV] is dedicated to study the sensitivities to 6*13 (in presence 
of NSI) and to and , which are achieved mainly through the Pg^ and Pgr channels; 
section |V] is devoted to the study of e^r and eaa, whose sensitivities are achieved mainly 
through the P^^ channel; in section I VI I we study the CP discovery potential of the three 
setups in presence of NSI; and finally, we conclude in section IVlTl In the Appendix we show 
approximate expressions for the oscillation probabilities Pg/^, Pgr, P^i^i and P^^r in matter with 
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constant density. 



II. NEUTRINO OSCILLATIONS WITH NON-STANDARD INTERACTIONS 

NSI can be studied from a top-down or a bottom-up approach. In the top-down approach, 
a given NP model is studied, and the corresponding set of low-energy effective operators are 
derived systematically. Being these operators derived from a fundamental theory, their 
coefficients are related and (usually) stringent bounds exist between them. In the bottom- 
up approach, on the other hand, all the effective four-fermion operators which can affect 
neutrino oscillations are included in the analysis, and the experimental bounds are used 
to constrain their coefficients independently of the model. This latter approach is clearly 
model independent, in the sense that no assumption is made regarding the model of NP 
behind and, consequently, the possible relations between operators. Bounds obtained in this 
way are looser than in the top-down approach, but apply to a wide variety of high-energy 
extensions of the Standard Model. 

In the bottom-up approach, NSI that modify neutrino production, propagation and de- 
tection processes must be included Such a large number of new parameters in 



the analysis, however, makes it extremely difficult to extract any useful information from the 
results. It is a standard strategy, thus, to separate the study of NSI in neutrino propagation 
in matter from the NSI effects in production and detection. The latter can be studied using 
near detectors [sl, |59|0. This is the strategy that we follow in this paper. We introduce all 



the operators that modify neutrino propagation in matter at once, adding nine new parame- 
ters (six moduli and three phases) to the existing i^SM parameter space. A complete analysis 
with such a huge number of parameters is extremely demanding, from both the numerical 
and analytical point of view. For this reason the effects of NSI in matter propagation have 
been widely explored in the literature turning on only one new parameter in the analysis, 
or two at most (one modulus and one phase) liH-Si. In this work, we try to achieve two 



main goals: (1) We attempt a first complete phenomenological analysis of the potential of 
HENF to constrain all the NSI parameters which can contribute to propagation in matter; 
(2) We illuminate complicated correlations between the effects of NSI and uSM CP violating 
phases. We believe that the complete treatment of effects of NSI in matter propagation is 
an important step toward a better understanding of possible physics beyond the uSM. 

In the model independent approach, constraints on NSI parameters in propagation are 
very mild, generically at O{10^^) or even order unity (giI . However, from the theoretical 



point of view, such large values of the NSI parameters are not really expected. This is easily 
understood if one tries to find a model of NP responsible for NSI effects without enlarging 
the low-energy SM particle content. Effects of NP at high energies manifest at low energies 
through an infinite tower of non-renormalizable effective operators of dimension d > A which 
are invariant under the SM gauge group. These are weighted by inverse powers of the NP 
scale A: 

= CsM + y 5/:'=' + ^ 6C''=^ + . . . , (1) 
A A^ 

where Csm is the SM Lagrangian which contains all SU{3)c x SU{2)l x f/y(l) invariant 
operators of dimension c? < 4. The factors 1/A.'^~^ appear to suppress the effective oper- 

^ Notice, however, that near detectors will necessarily put bounds on the combination of NSI in production 
and detection processes. Several near detectors with different sources and/or target materials could be 



used to disentangle them (see Ref. [60[) 
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ators which generate neutrino masses and produce NSI effects at low energy. Moreover, 
the necessary requirement of gauge invariance of the new operators under the SM gauge 



group [35|, |61|, |63H65| leads to another remarkable point: the effective operators which gen- 
erate neutrino NSI may be tightly related to their analogues in the charged lepton sector, 
which are much more constrained experimentally. Therefore, it is very hard to construct a 
feasible model giving large NSI effects in neutrino oscillations. 

Following the model independent approach, NSI in neutrino propagation (from here on, 
we will refer to them simply as NSI) are described through the inclusion of the following 
four fermion effective operators: 

<5£nsi = -2V2GFY.ei';{i7^YPLyp) Ql.Pf) , (2) 

f,P 

where Gi;' is the Fermi constant, / stands for the index running over fermion species in 
the Earth matter, / = e,u,d, P stands for the projection operators Pl = |(1 — 75) or 
Pji = |(1 + 75), and a, (3 = e,/i, r. Notice that from neutrino oscillations we have no 
information on the separate contribution of a given operator with coefficient e^^, but only 
on their sum over flavours and chirality. The effects of these operators appear in the neutrino 
evolution equation, in the flavour basifl as: 










^0 

















A21 





Vr ) 






^0 





A31 



1 + ee 



e, 



e* e* e 









)] 


1 








K-rJ 



(3) 



where = Am? /2i?, U is the lepton flavor mixing matrix, A = I^^Gpne and eap, = 

(l/-e) E/,pn,.S, with the /-type fermion number density. The three diagonal entries 
of the modifled matter potential are real parameters. Only two of them affect neutrino 
oscillations: we will consider the combinations eee — Ctt and e^^ — e^r, subtracting e^r x I from 
the Hamiltonian. The three complex NSI parameters ee^, e^T and e^r will be parametrized 

as|e„^ = |e„/3|e-'^-^. 

In order to understand the impact of different NSI parameters in various oscillation chan- 
nels it is useful to obtain approximate analytical expressions for the oscillation probabilities. 
In Ref. approximate formula, were derived for all the oscillation probabilities up to order 

[c" for the golden channel) by making a perturbative expansion in Am2i/Am|]^ = e and 
ea/3 ~ 6^13 ~ £. In the Appendix we present the approximate expressions for Pefii Peri PfifJ. 
and P^r up to second order in e expanding in Amli/Amg]^, 6^/3, 6^13 and 6623 = 6*23 — vr/4, 
too. 

Let us review very briefly the main conclusions which can be extracted from this analytical 
study: 

• Up to second order in e, Pe^ and Per depend only on and eer but not on the rest of 
the NSI parameters (see Eqs. lAll and IA2l in the Appendix). The precise determination 
of these two NSI elements is only possible in the golden and the silver channels. 
However, it is well-known that already in the uSM case a combination of data, either 
from different oscillation channels or different baselines, is needed in order to avoid the 
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If production or detection NSI were present, though, the effective production and detection flavour 
eigenstates would not coincide with the standard flavour ones 



66 



^ This is the prescription used in the MonteCUBES software. In the section devoted to CP violation, 
though, the prescription is precisely the opposite, e^^ = |eQ.^|e*'^°'^ . 
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well known degeneracy problem [68H70[. This problem is even more difficult to solve 
in presence of NSI, because four extra parameters (2 moduli and 2 phases) appear 
simultaneously in the golden and silver channels and severe correlations are expected 
to exist, not only between NSI parameters but also between them and the z/SM ones. 

• The P^f^ and P^,- oscillation probabilities show a leading 0{e) dependence on the real 
part of (which provides a very high sensitivity to this parameter), in addition to 
the usual quadratic dependence on and eer as in the golden and silver channels. 
On the other hand, the sensitivity to the imaginary part of e^^ is expected to be 
much worse, since it comes only through O(e^) terms in the probability. The depen- 
dence on the diagonal combination (e^^ — err) appears in P^^ and P^r at O^e"^), too. 
Terms proportional to 56'23(e^^ — e,-,-) lead to important correlations between these two 
parameters. 

• The dependece on (e^^— e^r) and e^^ in P^^ and P^^ is the same. Therefore, the z/^ — )■ Ur 
channel may be useful only because it adds further statistics at the detector. However, 
the sensitivities to eaa and e^r are not limited by statistics, since the disappearance 
channel alone already provides enough events at the detector. As a consequence, the 
sensitivities to these parameters are mainly achieved through the P^^ channel. 

• The dependence on the diagonal combination (eee — Crr) appears at third order in e in 
the oscillation probabihties. Therefore, it is hard to expect a good sensitivity to this 
parameter. Moreover, as we can see in Eq. ([3]), when all NSI parameters vanish except 
for the combination {eee — Ctt), A {eee — ^tt) can be interpreted as a small perturbation 
on the standard i/SM matter effect. Therefore, our sensitivity to {eee — ^tt) will be 
ultimately limited by uncertainties of the earth matter density. 

In view of the features listed above, in the following we are going to distinguish two 
different groups of oscillation parameters: (i) 6*13, ee^ and egr, that will be studied in Sec. HV] 
and (ii) (eee — e^r), (e/x^t — Crr) and e^r, that will be studied in Sec. |Vl This classification is the 
natural consequence of the fact that, in practice, Pe^ and Per are sensitive to (i), while the 
sensitivity to (ii) comes mainly from P^^. The only possible exception to this classification 
is (eee — ^tt)-, for which the golden channel also plays an important role. Because of this 
structure, strong correlations between (i) and (ii) are not expected, as it has indeed been 
found in our numerical simulations. For this reason, we will study in Sec. |Vl] the CP 
discovery potential of the HENF for all the parameters belonging to (i) simultaneously (as 
strong correlations are expected between them), whilst neglecting parameters belonging to 
(ii). 



III. THE STATISTICAL APPROACH AND THE SETUP 

In this section, we ffist introduce the statistical approaches used to perform the numerical 
analyses of Sees. HV] and |V] (see Sec. IIII Ap . and of Sec. |Vl] (see Sec. IIIIBp . After that, we 
recall the input values for the z/SM parameters and the current bounds on the NSI parameters 
(that are included as priors in the algorithm) and describe the marginalization procedure 
(Sec. IIII Cp . Eventually, we define the three HENF experimental setups which we are going 
to study (Sec. UTTDj) . 



6 



A. The statistical procedure used in Sees. IIVI and [V] 



It is well-known that, in order to sample a iV-dimensional parameter space through 
grids with n samplings per parameter, a total of (9(n^) evaluations of the expected number 
of events are required. When only three-family oscillations are considered, the computation 
can become heavy (if all i/SM parameters are taken into account) but is still affordable within 
the standard frequentist approach. When the NSI parameters are also taken into account, 
however, the number of parameters to be fitted simultaneously increases considerably and 
the computation time required to perform the standard minimization procedure becomes too 
large. A different approach must therefore be used if we want to sample a huge number of 
parameters with limited computational resources. The way out is suggested by noticing that 
most of the points belonging to the grids that are computed in the standard approach are 
useless, as they are very far from the minimum. For this reason the standard technique 
used to sample large multi-dimensional manifolds is to rely on efficient (either deterministic 
or stochastic) algorithms that search for the global minimum and then start to sample 
the region near the minimum to determine its size and shape. Most of the algorithms 
used fall into the category of Markov Chain Monte Carlo (MCMC): using these class of 
algorithms, the number of evaluations required for the algorithm to converge and sample 
properly the desired distribution grows polynomially with A^, 0{N^), with k some integer. 
We have followed this approach to scan the NSI parameter space in Sees. IIVI and IVl using 
the MonteCUBES ("Monte Carlo Utility Based Experiment Simulator") software |47| that 
contains a C library plug-in to implement MCMC sampling into the GLoBES Tlf package. 



It, thus, benefits from the flexibility of GLoBES in defining different experiments while 
implementing an efficient scanning of large parameter spaces. 

Parameter determination through MCMC methods are based on Bayesian inference. The 
aim is to determine the probability distribution function of the different model parameters 
6 given some data set d, i.e., the posterior probability P{6 \ d). From Baye's theorem we 
have: 

The likelihood L^iO) = P{d \ 6) is the probability of observing the data set d given certain 
values of the parameters 6. The prior 71(6) = P{6) is the probability that the parameters 
assume the value 9 regardless of the data d, that is, our previously assumed knowledge of 
the parameters. Finally, the marginal probability M is the probability P{d) of measuring 
the values d. It does not depend on the parameters 6, and therefore it can be regarded 
as a normalization constant|j. Notice that the functions defined in GLoBES provide 
the logarithm of the likelihood of the data d following a Poisson distribution normalized 
to the distribution with mean d. Therefore, the actual probability density sampled by 
MonteCUBES is the posterior probability P{6 \ d): 



V = exp 



exp 



where Xp(^) = — 21n7r(^). This probability distribution is equivalent to a Boltzmann weight 
with temperature T = 1 and energy E = x^i^) + 



^ In the limit of infinite statistics, it can be shown that the Bayesian probability distribution is maximized 
by the same set of parameters 9 that minimize the function in the frequentist approach. 
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We have used ten MCMC chains in all our simulations. The convergence of the whole 
sample improves as R ^ 1, with R being the ratio between the variance in the complete 
sample and the variance for each chain. We have checked that the chains have reached 
proper convergence in all cases better than R — 1 = 2.5 x 10~^. 

A typical problem when a minimization algorithm different from the complete computa- 
tion of the mult i- dimensional grid is applied is the possible presence of local minima or of 
multiple global minima ("degeneracies"). In both cases, if the minima are deep enough the 
algorithm will get stuck there and sample a region that does not correspond to the global 
minimum or will not be able to identify the presence of degenerate minima. The Monte- 
CUBES package includes a method to identify local minima by increasing the temperature 
T of the chain so that the likelihood is modified to "P oc V^^^ . This procedure flattens the 
likelihood distribution, making it possible for the chains to jump from a local minimum to 
another. The temperature and step sizes are then decreased in successive steps and thus the 
different chains get stuck around different minima, unable to move through the disfavored 
regions when T is too low. After this, the points where the different chains have stopped 
are compared to decide how many different minima the chains have fallen into. Finally, new 
steps are added with the correct length in the direction between the degeneracies. As a 
result, the algorithm is able to jump between minima and sample all of them properly. 

For the implementation of the NSI probabilities in matter, we use the non-Standard 
Interaction Event Generator Engine (nSIEGE) distributed along with the MonteCUBES 
package. 

The definition of confidence level (CL) in a multi-dimensional MCMC algorithm need to 
be clarified, as it approaches the standard definition only in the limit of infinite statistics. 
What is done in practice with a MCMC is to generate a given number of points distributed 
stochastically in the multi-dimensional parameter space around the global minimum (min- 
ima). After that, projections of that multi-dimensional "cloud" of points are performed 
onto any desired plane chosen accordingly to the variables under study in the analysis. The 
projection over a given plane corresponds to marginalization over the parameters that are 
no longer retained. After projection, the two-dimensional plane is divided into cells and the 
number of points falling into each cell is computed. For small enough cells, the resulting 
two-dimensional histogram can be approximated by a smooth surface, for which slices can 
be drawn for the desired CL (68%, 90%, or 95%, in this paper). 

B. The statistical procedure used in Sec. IVII 

To explore the CP discovery potential we need a different statistical procedure from the 
one defined in Sec. lIIl7\l The reason is the following: a MCMC, as described above, explores 
the region which is close to the global minimum (or to degenerate minima) sampling with 
good accuracy the distribution around that point (s). This is the right procedure to follow 
if we are exploring the sensitivity that a facility has to some particular observable. Consider 
the particular case of the sensitivity to ^13 in the (^13,5) plane with marginalization over 
the rest of z/SM and NSI parameters (see Sec. lIVp . In this case, the MCMC algorithm scans 
the multi-dimensional surface corresponding to a given choice of the input parameters (with 
the particular choice ^13 = for 6*13) and a contour at a given CL of the region compatible 
with vanishing 61^ is drawn, after projecting over the (^13,5) plane. When we compute a 
discovery potential, on the other hand, we first fix the parameters to be tested and draw the 
corresponding CL contours. Then, we check if the condition we want to fulfill is satisfied 



8 



or not at a given CL (in the case of the CP discovery potential, the condition is that the 
contours drawn for a given set of CP violating input parameters do not touch any CP 
conserving point of the parameter space). Eventually, we repeat the procedure again and 
again varying the input parameters. If the grid density is large enough, the distribution of 
the input parameters that satisfy the required condition is smooth and a"CL contour" can 
be drawn. 

If we were to use the MCMC approach to compute a discovery potential, then, we should 
run the algorithm as many times as the points in the grid that we want to test. In this 
case, the total time required to compute the discovery potential goes as n^^ x A^^, with 
n the number of points to be tested for one parameter and A''^ the dimensionality of the 
grid. If n is large the MCMC cannot be used and the standard frequentist approach must 
be adopted instead. The drawback of the frequentist approach is that, in order to keep 
the computational time from being rapidly divergent, we cannot marginalize over the whole 
i/SM and NSI parameter space. For this reason, in Sec. |Vl] we will not marginalize over 
atmospheric and solar parameters and will consider fixed inputs for the NSI parameters. 

The procedure that has been used in this work to determine the CP-discovery potential 
is outlined below: 

1. We first compute the number of expected events at the detector(s): A^(6'i3, {0})? where 

2. After having computed grids of number of events as a function of (^13; {0}), we com- 
pute the as follows: 

polarities, bins 

with / an overall systematic error. In all the plots given in this section, we assume 
the overall systematic error for the MIND detector as = 0.02, and the one for the 
r-signal as fr = 0.05. No background has been considered to compute Eq. ([5]). 

3. We compute, then, for any input (^13; {0}), the function defined in Eq. at the 
eight CP-conserving (CPC) points: 

{0}cpc = (0, 0, 0); (0, 0, vr); (0, vr, 0); (vr, 0, 0); (0, vr, vr); (vr, 0, vr); (vr, vr, 0); (tt, tt, tt) , 

taking the smallest value found. Using this procedure, we obtain the five- 
dimensional surface: 

XcPciOis, 0,s; {0}) = min {x'iO^s, ^13; {<P}cpc, {0})) (6) 

iViCPC 

This procedure generalizes to the case of three simultaneously active phases the pro- 
cedure outlined in Ref. where only S and (j)eT were considered. 



4. We distinguish, then, between two cases depending on the value of 6 



13- 



The first possibility is that ^13 is already measured by the time the HENF is 
built. In this case, we can use the function computed at ^13 = ^13 to see 
the region of the phase parameter space where CP violation can be distinguished 
from the CP conservation hypothesis. 
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The second possibility stands for a very small (or even vanishing) ^13. In this case, 
it is also necessary to marginalize over ^13 since possible CP-conserving solutions 
can be found for a given CP-violating input (^13; {0}) at a different 6^13 (what in 
the standard three-family oscillation scenario is called an "intrinsic degeneracy" 



68j). In order to take these degeneracies into account, we minimize the over 



^'13: 

xIcpM^; {0}) = min (xcpc(^i3, ^13; {0})) • 
013 

However, notice that in this case the only information we have on i9i3 is an 
upper bound, ^13 < 3°, approximately. Therefore, marginalization over ^13 in the 
allowed range is also required here: 

xle,cpc({^}) = (Xe,c7Pc(^i3; {0})) • (7) 

5. Eventually, we draw the three-dimensional surfaces corresponding to = 11.34. 
These contours represent the area of the phases parameter space in which CP vio- 
lation can be distinguished from CP conservation at the 99% CL for 3 d.o.f.'s. Results 
will be shown for both cases in which ^13 is known, using Eq. or unknown, using 
Eq. dZD. 



C. Input parameters and marginalization procedure 

Unless otherwise stated, the input values taken in this work for the atmospheric and 
the solar parameters are: Amgi = 7.59 x 10"^eV^, 612 = 34°, Amg^ = 2.45 x 10"^eV^, 
623 = 45.5° jloj. In all the simulations, the matter density has been taken according to the 
PREM density profile [72[ assuming a 5% error. The sign of the atmospheric mass difference, 
sgn{Amli), has been chosen to be positive throughout the paper and marginalization over 
it will not be considered. 

In all the simulations presented in Sees. |IV]and|V]we have marginalized over the whole set 
of z/SM parameters. A gaussian prior distribution centered on the input values given above 
with variance a = 0.08(0.03) has been assumed for the atmospheric (solar) parameters. On 
the other hand, marginalization over ^13 and 6 has also been performed assuming a flat prior 
distribution. 

In Sees. IIVI and |V] we have also marginalized over the NSI parameters eai3, that is, 
over both their moduli {eapl and phases 0q,/3. Gaussian priors, in agreement with the bounds 
computed in Ref . 6^ , are taken into account for all the moduli of the NSI parameters around 
their input values, which have been set to zero throughout the next two sections. Notice that 
for the NSI phases no prior knowledge has been taken into account {i.e., 7r{6) = 1), since 
we do not have any information about these phases yet. We will refer the above procedure 
involving uSM and NSI parameters as the "standard marginalization procedure" hereafter in 
this paper. 



D. Neutrino Factory setups 

At a NF, intense z/g and z/^ beams are available as decay products of muons (with both 
polarities) circulating in storage ring(s). As a consequence, a total of twelve different oscilla- 
tion channels could in principle be studied. The International Design Study for a Neutrino 
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Factory |50[, as already mentioned, has undertaken the task of defining the optimal setup 
to have good sensitivity to ^13, S, and to the neutrino mass hierarchy, sign(Am23). This 
resulted in what we will refer to as the IDS25 setup hereafter: a HENF with a muon beam 
energy of 25 GeV, and 10^^ useful muon decays per year aimed at two identical 50 kton 
MIND detectors located at two different baselines to look for appearance event^. The 
detector located at L = 4000 km (which will be referred to as "intermediate baseline" from 
now on, to distinguish it from the short baseline where a near detector will be locatec^ ) is 
optimized to have sensitivity to the CP-violating phase S. It, however, suffers from a severe 
degeneracy problem 6^70|, that can be solved locating an additional "far" detector at the 



so-called "magic" baseline (7500 km) [68|,|73|. This second detector increases significantly 



the potential of the NF to measure the mass hierarchy, taking advantage of matter effects. 

If a HENF is to be optimized to detect effects of NSI, several issues must be understood: 
which energy and baseline would be the best; how and to what extent the synergy between 
two detectors help; and, if there are any ways of optimization in order to achieve good 
sensitivities to both z/SM and NSI parameters. Some of these problems were addressed in 
(4o| . concluding that a setting similar to IDS25 but with higher muon energies (such as 50 
GeV) would be preferred to look for NSI. This is easily explained by the fact that, since NSI 
in propagation are introduced as an effective matter potential, an increase in the average 
neutrino energy will improve the relative significance of the NSI with respect to the leading 



standard oscillations in vacuum. This was indeed confirmed in Ref. j4l| where, however. 



it was found that the improvement with respect to the 25 GeV setup was not very large 
(see, also, Ref. j42|). All of these works were performed within the ad-hoc assumption 
of having only one e^^ at a time. One of the goals of this paper is indeed to check if 
the results obtained in these analyses survive when correlations between the various NSI 
parameters are taken into account. Armed with MonteCUBES, we examine this problem 
by comparing the sensitivity to NSI of the IDS25 and of a variant of the IDS setup with the 
same detectors but with parent muon energy E^^ = 50 GeV when all NSI parameters are 
turned on simultaneously. The new setup, for obvious reasons, will be called as IDS50. 

We will also check the potential of the two detectors to reduce the strong correlations 
between z/SM and NSI parameters and within different NSI parameters, to test if this detec- 
tor combination can be optimized or not for NP searches. For this reason, the performance 
of the two setups above, in which two identical detectors looking for z/e — ?■ z/^ and z/^ — ?■ z/^ 
oscillations are located at different baselines, will be compared with the performance of a 
HENF setup in which two different detectors, one of which is equipped with r- identification 
capability (that could, thus, profit of the z/g — )■ and z/^ — )■ Ur channels), are located at 
one single baseline at 4000 km. Initially, L = 2000 and L = 3000 km were also considered 
as alternative baselines. They give similar results, although the 4000 km performs slightly 
better, and therefore will not be considered here. This setup will be called 1B50. Notice 
that the parent muon energy for the 1B50 setup has to be large enough to overcome the 



smallness of the u^-N cross-section due to the r production threshold below 4 GeV [7^ and 
get larger statistics at the detector. For this reason, also for this setup we fix = 50 
GeV. Notice that the advantage of aiming at one site is two-fold: on one side, we avoid the 
technical difficulties of aiming one of the beams at a detector located al L = 7500 km (with 
a tilt angle of the storage ring of ~ 36° [i^); on the other side, as only one storage ring is 



needed all muon decays are aimed to the same site, therefore doubling the statistics at the 



^ Notice that, in the latest IDS design, the MIND located at the intermediate baseline is 100 kton. 
® We will not consider such a near detector in our analysis, though. 
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detector. 

The characteristic features of the three setups (IDS25, IDS50 and 1B50) are resumed in 
Tab. [U where we remind the parent muon energy, the detectors location and technologies 
and the neutrino flux aiming at each detector per year. For all setups we consider 5 years 
of data taking for each muon polaritjfl. 





IDS25 


IDS50 


1B50 




25 GeV 


50 GeV 


50 GeV 


Di 


MIND@4000 km 


MIND@4000 km 


MIND@4000 km 


D2 


MIND@7500 km 


MIND@7500 km 


MECC@4000 km 


$1 


5 X 10^° 


5 X 10^° 


1 X 10^^ 


$2 


5 X 10^° 


5 X 10^° 


1 X 10^^ 



TABLE I. The characteristics of the three considered setups. Prom top to bottom: the parent 
muon energy E^; the technology and location of the two detectors, Di and D2; the number of 
useful muon decays per year aiming at each of the two detectors, and <I>2- 



The characteristic features of the two types of detector are summarized in Tab. HTl For 
technical details on these parameters, we address the interested reader to Refs. j48[ and 



451,178 





a{E) 


fs 


fB 


Mass 


MIND 


0.55/B 


2.5% 


20% 


50 kton 


MECC 


0.2^ 


15% 


20% 


4 kton 



TABLE II. Main characteristics of the two detectors technologies. From left to right: energy 
resolution, cr{E); systematic error over the signal, fs; systematic error over the background, /s; 
detector mass. 

In the analyses of Sees. [IV] and |Vl data have been distributed in bins of the reconstructed 
neutrino energy with the following size: AEi, = 1 GeV for G [1, 10] GeV; AEi, = 2.5 
GeV for E^ e [10, 15] GeV; AE^ = 5 GeV for E^ > 15 GeV. This binning applies to aU 
setups and detector technologies. In Sec. IVTl on the other hand, data have been distributed 
in bins of equal size: AE^ = 5 GeV for the IDS25; AE^ = 10 GeV for the IDS50, 1B50 
(both for the MIND and the MECC technologies). 

^ We assume that the experiment is run on the two polarities separately. This means that we are considering 
a total number of useful muon decays per baseline and polarity of 5 x 5 x 10^° = 2.5 x 10^^. Notice that 
this is completely equivalent to consider 10 years of data taking per polarity but with 2.5 x 10^'' useful 
muon decays per baseline, year and polarity. This last option corresponds to the setup where muons of 
both polarities are circulating at the same time in the decay ring(s). 
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The efficiencies for the MIND and the MECC detector technologies as a function of the 
reconstructed neutrino energy are shown in Fig. [H The i/^ identification efficiency at MIND 



has been taken from Ref. [48||. The Ur identification efficiency at MECC corresponds to the 
efficiency of the ECC for the silver channel i/g as computed in Ref. [75| multiplied by 

a factor of five to take into account the capability of the MECC to look for taus not only 
through their decay into muons (as for the ECC) but also into electrons and hadrons (see 



from Refs. 76, 77 



a detailed discussion in Ref. ||45| regarding this point). The cross-sections have been taken 
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FIG. 1. Efficiency of the MIND (blue triangles) and MECC (red circles) detectors as a function 
of the neutrino energy. The MIND efficiency has been taken from Ref. 48|]. The MECC efficiency 
corresponds to the ECC efficiency [75| multiplied by a factor of five (see text for details). 



Table Uni shows the number of events in the golden (silver) channel per ktonxyear at a 
MIND (MECC) detector with perfect efficiency, located at 4000 km from the source, for a 
50 GeV NF. Normal hierarchy has been assumed, and results are shown for b = ±90°, and 
for two different values of ^13 = 0, 3°. In order to illustrate the effect when NSI are included 
in the analysis, we show the total number of events at the detector also in presence of NSI, 
The two NSI CP violating phases have been set to zero. As it can be seen 



-efi 
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from the table, the number of events in the silver channel is very small for vanishing 9 



13- 



IV. SENSITIVITIES ACHIEVED MAINLY THROUGH THE i^f, AND Vr 
CHANNELS 

As we stressed in Sec. [Tll the sensitivity to NSI parameters comes from different oscillation 
channels depending on the considered parameter. In this section, we study the sensitivities 
to and eer which would be achieved mostly through the golden (and, to a lesser extent, 
also through the silver) channel for the three setups under study. Since the sensitivity 
to the z/SM parameter ^13 (which is the key to the measurement of 6 and of the mass 
hierarchy, too) is also achieved through the same oscillation channels, we will examine first 
the question of how and to what extent the inclusion of NSI affects the sensitivity to 
(Sec. IIV A|) . We will, then, study the sensitivities to the moduli of ee^and eer as a function 
of their respective CP violating phases, through which some features of the synergy between 
two detectors/baselines will be illuminated depending upon the settings (Sec. IIV Bp . To 
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6 


Channel 


^13 = 
ea/3 = 


^13 = 
— £er — 10 


^13 = 3° 
Cq/? = 


^13 = 3° 

^e^l — £er — 10 


+90° 


T^ti 


2.75 


44.93 


66.70 


108.54 


Ve Vr 


1.08 


7.97 


15.63 


22.60 


90° 


Ve l^fi 


2.75 


44.93 


41.69 


83.96 


Ve Vt 


1.08 


7.97 


23.77 


30.56 



TABLE III. Total number of events per year for the golden (silver) channel, measured at a 1 kton 
MIND (MECC) detector with perfect efficiency located at L = 4000 km from the source, for a 50 
GeV NF. Normal hierarchy has been assumed. Results are presented for 6 = ±90° and for two 
different values of ^13 = 0,3°, with and without including NSI effects in the golden sector. The 
two NSI CP violating phases have been set to zero: (pe^i = (per = 0°. 

show the effect of the correlations between the NSI parameters and with ^13, the standard 
marginalization procedure defined in Sec. IIII CI will not be always used. We will specify 
in each case the procedure adopted. We remind the readers that sgn{Aml^) is kept fixed 
throughout this paper. This point has to be kept in mind when interpreting the results. 

To conclude the preamble, we describe the layout of the figures: in figures with three 
columns, the left, middle and right panels correspond to the results for the IDS25, IDS50 
and 1B50 settings, respectively. Red, green and blue lines correspond to 68%, 90% and 95% 
2 d.o.f.'s CL contours, respectively. Whenever we depart from the standard format we will 
give a note in the caption of the corresponding figure to specify the layout. 



A. Impact of the NSI on the measurement of ^13 

In Fig. [2] we show the sensitivity to 613 as a function of 6 when the NSI parameters, 
in addition to the z/SM ones, are also taken into account during marginalization. In top 
panels, marginalization over eaa = e, /i, r), ee^and eg^is performed; in bottom panels, the 
marginalization procedure is done over eaa and e^^-. As a reference, we also present the 68% 
CL sensitivity to 6*13 obtained without considering NSI in the analysis, represented by the 
dotted black line^. It can be clearly seen that the impact of the presence of NSI degrees of 
freedom on sensitivity to ^13 is much more significant in top panels than in bottom panels. 
When marginalization is performed over ee^and eer, the degree of sensitivity loss ranges from 
a factor 3 (IDS50) to almost an order of magnitude (1B50) with respect to the uSM result. 
On the other hand, the effects of marginalization over eaa and e^^ are quite mild, leading to 
a sensitivity loss of a factor of 3 (1B50), at most. We have checked that the above sensitivity 
loss comes from the marginalization over eaa, and not over e^^ which is effectively decoupled 
from 6'i3. Approximate decoupling between the two parameter sets {613, e^^ , ) and {eaa-, 



We have checked that our results for the i^SM 0i3-sensitivity are in reasonable agreement with those 



reported in the literature (see, for example, Ref. 17j ) 
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e^r) is consistent with the expectation from the perturbative analysis (67[ (see Appendix 
as explained in Sec. [Ill in a nutshell, the golden (and silver) channel oscillation probabilities 
(that dominate the sensitivity to 613) only depend on ee^^and eerUp to second order in e. 




IQ-'- 10-' lO-* 10"' 10-5 10-4 jo-6 10-5 10-4 

sin2(20,3) sin\2e,3) sin2(2ei3) 



FIG. 2. 68%, 90% and 95% CL contours for the sensitivity to ^13 as a function of 5 for the case with 
NSI, compared to the 68% CL contour for the sensitivity to ^13 in the absence of NSI (represented 
by the black dotted line). Marginalization was performed over the vSM. parameters, the diagonal 
NSI parameters eaa and either Cg^i and egr (top panels) or e^,- (bottom panels). The left, middle, 
and right panels show the results obtained for IDS25, IDS50, and 1B50 setups, respectively. 



A careful comparison between the three upper panels reveals an interesting feature: the 
severe impact on the sensitivity to ^^13 observed at the 1B50 setup is largely (moderately) 
overcome in the IDS50 (IDS25) setups. Sensitivity to 6is at the two baseline settings is 
robust against inclusion of NSI because they probe generalized matter effects at two different 
distances, a particular type of the synergy between the intermediate and far detectors |40 . 



In summary, in spite of the fact that the 1B50 setup apparently yields the best sensitivity 
to 6*13 in absence of NSI (which is likely to be due to the doubled flux with respect to the 
IDS25 and IDS50 setups), its worsening after marginalization over and eer , strongly 
correlated with 6*13, is much more severe due to the lack of the magic baseline detector. 



B. Sensitivity to and 

In Fig. [3] (jl]) we present the sensitivity to |ee^| (|eer|) for ^13 = (upper panels) and 
= 3°; 6 = —7i/2 (lower panels) as a function of 0e/x {(per), respectively. For all panels, 
the standard marginalization procedure defined in Sec. IIII CI is carried out, albeit neglect- 
ing marginahzation over e^,- (which is totally uncorrected from |ee^| and |eer|, as we have 
checked) . 

We first discuss the results in Fig. [31 The most important remark is that both high- 
energy setups (IDS50 and 1B50) present similar performances, with sensitivities that are 
much better than that of IDS25 both for ^13 = and ^13 = 3°. The fact that both IDS50 
and 1B50 give sensitivities to \eefj,\ that are extremely similar implies that the improvement 
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FIG. 3. 68%, 90% and 95% CL contours for the sensitivity to as a function of (fieij. for ^13 = 
(upper panels) and ^13 = 3°; (5 = — 7r/2 (lower panels). Marginalization has been performed over 
the z^SM parameters, the diagonal NSI parameters eaa and eer- The left, middle, and right panels 
show the results obtained for IDS25, IDS50, and 1B50 setups, respectively. 



with respect to the IDS25 is due to the increase in energy, in agreement with our expectation 
for preferring higher energies for NSI searches. It can also be seen that the dependence of 
the sensitivity to \eep\ on ip^^ is stronger for ^13 7^ 0° (bottom panels) than in the case of 

^13 = 0. 

In Fig. m we can see that the results obtained for the sensitivity to |eer| are quite dif- 
ferent from those found for |ee^|. For vanishing ^13, it can be seen that the IDS50 setup 
is better than the IDS25 and the 1B50 by a factor of ~1.7 and a few, respectively. When 
^13 is increased to ^13 = 3°, the sensitivity to eer becomes worse by a factor of ~2 or so, 
independently on the setups: the relative performance of the three setups remains almost 
the same as in the case of ^13 = 0. Notice that for the sensitivity to |eer|, unlike for \ee^\i 



the synergy between the two detectors plays a key role |40H42| . improving the sensitivity up 
to an order of magnitude when a second detector at the magic baseline is considered. An 
additional improvement is achieved due to the increase in energy, as expected. 

The oscillation probabilities presented in Appendix |X] can help us to understand further 
the results presented in Figs. Eland ID In the first line of P^^ in Eq. (lAip we can observe that 
the two terms proportional to ee^i appear with the same sign, while the ones proportional 
to eer have opposite sign and tend to cancel. We have checked that, in the energy range 
relevant for the three setups under consideration, the coefficient of the ee^ term can be one 
order of magnitude larger than the corresponding coefficient of the eer term. For this reason, 
the golden channel is more sensitive to ee^than eer, as it can be seen by comparing Figs. [3] 
andm For the same reason, the eg^ sensitivity is only mildly affected by marginalization over 
eer, whereas the sensitivity to egris strongly affected by the marginalization over ee^. This 
explains why the detector at the magic baseline plays an important role for the sensitivity 
to I eer I, whereas for ee^ the energy is the key parameter independently of the number of 
baselines. Notice that the features of Per are quite the opposite: terms proportional to 
eer add up, while those proportional to ee^^ tend to cancel (see the first line of Eq. flA2|) in 
the Appendix). One could naively think, then, that a better sensitivity is expected for eer 
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FIG. 4. 68%, 90% and 95% CL contours for the sensitivity to Cer as a function of c^er for ^13 = 
(upper panels) and ^13 = 3°; 6 = —tt/2 (lower panels). Marginalization has been performed over 
the uSM parameters, the diagonal NSI parameters eaa and ee/^. The left, middle, and right panels 
show the results obtained for IDS25, IDS50, and 1B50 setups, respectively. Notice the different 
scale for the lower right panel, for which the sensitivity is much worse than for the rest of setups. 



instead of eg^in the 1B50 setup due to its ability to detect r's. This is not the case, however, 
because the statistics in MECC is poor compared to that in MIND. The good sensitivity 
observed at the 1B50 depends on the doubled flux at the MIND section of the detectoiQ. 

A last interesting remark can be drawn in the phase dependence of the results of Figs. [3] 
and m If we compare the sensitivity contours for ^13 = (upper panels) and 613 = 3° (lower 
panels), we can see a shift of locations of the sensitivity minima. This feature is a result of 
the complicated correlations between 6, (p^fi and (per in the golden channel probability. The 
key factor for this effect to take place is the CP-violating value we have chosen for S, which 
maximizes the effect. 



V. SENSITIVITIES ACHIEVED MAINLY THROUGH THE i/^ ^ i^^ CHANNEL 

As we have mentioned in Sec. [ITl sensitivity to the rest of the NSI parameters, i.e. the 
diagonal elements and e^.^ , mainly comes from the z/^ — z/^ channel^. In this section, 
we will study first the expected sensitivity to \e^r \ as a function of its CP- violating phase 0^.^ 
(Sec. IV Ap . The sensitivity to the diagonal NSI parameters will be studied next (Sec. IV BO . 
taking particular care to unveil the correlations between eaa, diz and ^23- 

A. Sensitivity to e^r 

In Fig. the sensitivity to e^,- for ^13 = is shown only for the IDS50 setup, since 
we have found remarkably similar results for the rest of setups under study. The standard 

^ We will see in Sec. I VII that the silver channel can, however, play an important role in the discovery of 
CP violation due to NSI. 

The — > Vt channel only increases the statistics at the detector. However, the sensitivity to eaa and e^r 
is not limited by statistics thanks to the disappearance channel. As a consequence, the v^j_ — > Vr channel 
is not very useful in this context since it does not -add any additional information. 



marginalization procedure is employed as usual. The most significant feature in this figure is 
the extremely high sensitivity to the real part of e^r (better than 10~^). The high sensitivity is 
driven by the leading NSI correction to the disappearence oscillation probability in Eq. (\A3\i : 

P^^ = P^l - M coscP^riAL) sin (A31L) + 0(8^) + ... (8) 

When e^T- is mostly imaginary, (p^^- ~ ±90°, we observe a significant sensitivity loss of more 
than an order of magnitude, as expected from the fact that the leading dependence on 
Im(e^T-) appears at O(e^) in the probability, as seen in Eq. (]A3p . We have explicitly verified 
that fixing eefi and Cgr during marginalization does not produce any appreciable change in 
the sensitivity to e^,-. This confirms the numerical results presented in Sec. IIV Bl where no 
correlation between (or eer ) and e^r was found. Eventually, we have checked that the e^r 
sensitivity does not vary significantly for nonzero ^13, in agreement with the discussion in 
Sec. HTOl 
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FIG. 5. 68%, 90% and 95% CL contours for the sensitivity to e^,- as a function of </>^t- for ^13 = 0. 
Marginalization has been performed over the z^SM parameters, the matter density and the rest of 
NSI parameters. The 50 GeV IDS setup has been assumed. 




B. Sensitivity to the diagonal NSI parameters 

In Fig. E] we show the sensitivities to the NSI diagonal parameters obtained with the 
IDS50 setup for different input values of ^13 and ^23- As for the e^^ case, the results are 
extremely similar for the other two setups under study, and hence their results are not 
shown. Top panels correspond to ^13 = 0, the bottom ones to ^13 = 3°. Left panels are 
obtained for 623 = 45°, in which the red, green and the blue lines stand for the 68%, 90% 
and 95% CL contours, respectively. In the right panels only the 95% CL contours are drawn, 
for 023 = 43° by the purple dashed lines and 623 = 47° by the black solid lines. The standard 
marginalization is adopted here too. 

First of all, we see that all panels show that the sensitivity to (cge — ^tt) is about an 
order of magnitude worse than the sensitivity to (e^^ — err)- This behaviour is in agreement 
with the fact that the leading dependence on the latter combination appears at 0{e'^) in 



the oscillation probabilities, while for the former one it appears at 0{e^) (see 67|). Notice, 



nonetheless, that the sensitivity to (eee- Crr) (approximately ~ 10% and ~ 20% at 95%CL, 
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FIG. 6. Sensitivity to (egg -Crr ) and (e^^ -Crr ) at the IDS50 setup. Top panels: ^13 = 0; bottom 
panels: ^13 = 3°. Left panels: 68%, 90% and 95% CL contours for §23 = 45°; right panels: 95% 
CL contour for 6*23 = 43° (purple dotted line) and 623 = 47° (black solid line) . Marginalization has 
been performed over 823, Am^j^, the matter density, ^13 and 5. 



2 d.o.f.'s, for ^13 = 3° and 0, respectively) is better than the sensitivity achieved at any 
other facihties considered in the hterature. It improves as ^13 increases, probably due to the 
effect of the golden channel, as we can see by comparing top and bottom panels. On the 
other hand we have observed, in agreement with the discussion in Sec. HIl that the (eee - Crr 
) sensitivity is mainly limited by the matter uncertainty, which has been set to 5% in our 
simulations. In other words, unless the PREM error on the matter density is improved, a 
~ 10% sensitivity to (eee - Ctt ) would be the limiting accuracy that could be reachec0- 

The impact of a non-maximal atmospheric mixing angle can be seen in the right panels 
in Fig. ini two narrow strips appear at both sides of the central region. By looking into the 
disappearance probability, P^^ in Eq. (]A3p . it is easy to realise that the sensitivity in the 
central region of the plots is driven through the term proportional to (e^^ — e,-^)^ in the 
disappearence channel, being the only one which does not vanish to order for maximal 
mixing in the atmospheric sector. The narrow bands at both sides appear as a consequence 
of the non-maximal input for the atmospheric mixing angle, and they are driven by the 
terms proportional to S923{e^^ — err) in P^in- Indeed, the fact that these two strips appear 
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An alternative method to constrain e^e, which is free from this problem and is complementary to our 



method, is to use solar neutrinos, whose sensitivity to Cee appears to reach ^ 20% at 1 aCL (1 d.o.f.) 79[. 
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both for ^23 = 43° and ^23 = 47° indicates the existence of an "octant" degeneracy between 
56'23 and (e^^- err)- 

It is also remarkable that no significant correlations among eaa, ^efi and eer have been 
found, as expected from the results presented in Sec. IIVBI 



VI. DISCOVERY POTENTIAL FOR CP VIOLATION IN THE {(pe^^fpeT,^) SPACE 

One of the most interesting aspects of any system involving NSI is the possible existence 
of multiple sources of CP violation. It is important to understand characteristic features 
of CP violation such as correlations between the phases or possible degeneracies arising 
between them. In this paper we focus on the study of CP violation associated with the two 
NSI CP- violating phases 0e/iand together with the standard z/SM phase 6, which appear 
in the golden and silver channels. Depending on the values of 6*13, |ee^| and |eer|, strong and 
complicated correlations are expected to take place between these three phases. The third 
NSI phase, (p^r, only appears in P^^ and P^r (see Appendix and is uncorrelated to the 
rest of the CP-phases. For this reason effects of CP violation in the fi — t sector will not be 
studied in this paper. 

For any realistic NP model giving rise to NSI at low energies, the effects are generally not 
expected to be larger than (9(10~^). For this reason, we will focus on "reasonable" values 
for |ee^| and leerl, in the range \eai3\ G [10^'^, 10"^], for which the correlations with 613 can 
still be large ^n. For NSI moduli smaller than 10~^, the effect of S always dominates and it is 
very difficult to detect CP violation due to NSI. On the other hand, we will see that a very 
interesting structure arises in the discovery potential in the three-dimensional parameter 
space for the input values we have considered for the NSI moduh. In particular, we have 
studied three cases: (a) both moduh are "small", |ee^| = \eer\ = 10~^, (b) both of them are 
"large", 10~^, and finally (c) an "asymmetric" case where \eefj,\ = 10~^; \eer \ = 10~^. Within 
these three cases, the last one is particularly interesting. Let us remind that, in the golden 
channel, plays a leading role while eer is subdominant (Sec. lIVp . Therefore, the difference 
by a factor of ten in their order of magnitudes triggers interesting three-fold correlations 
between 6, (pe^i and (per- 



A. Non-standard CP violation in the absence of i^SM CP violation 

The first question we address is whether it is possible to detect a new CP-violating signal 
due to NSI in the absence of standard CP violation. We study, therefore, the CP discovery 
potential (defined in Sec. IIII B|) in the case where the input value for the standard CP- 



violating phase, 6, is set to zero or tt. The two cases defined in Sec. IIII Bl are considered. 
The first possibility stands for a relatively large value of ^13, 613 = 3°, which would be already 
measured by the time the HENF is built. In this case, we can safely use Eq. to study 
the CP-discovery potential. The second possibility arises when no signal for a non-vanishing 
di3 is found by the ongoing and soon-coming neutrino oscillation experiments. In this case, 
we can only conclude that 613 < 3°, and we have to use Eq. ([7]) instead. These two cases 



The range of values considered for the NSI moduh in this paper differ significantly from that in Ref. [43| , 
where values for {eerl as large as unity were considered. In addition, the number of NSI parameters 
in the analysis is different as well. As a consequence, the comparison of the results obtained is not 



straightforward. Some qualitative features of the results obtained in Ref. [43| have been recovered, though 
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yield the best and the worst resultfl for the CP discovery potential in the {(pefj,, (per) space 
for 5 = and 180°. 
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FIG. 7. The 99% CL (3 d.o.f.'s) CP discovery potential in the {(p 

e flipper) plane for 5 — 0. From 

top to bottom: (|ee;,|, |ee^|) = (10"^, 10"^), (IQ-^, IQ-^) and (IQ-^, IQ-^). The yellow (light gray) 
regions have been obtained for ^13 = ^13 = 3°, Eq. ([6|). The cyan (dark gray) regions have been 
obtained after searching for intrinsic degeneracies in ^13 and then marginalizing over ^13, Eq. ([7]). 



We show in Figs. [7] and [8] the CP discovery potential in the {(pefj,, (per) plane for 5 = 
and 180°, respectively. From left to right we show results for the IDS25, IDS50 and 1B50 
setups, respectively. From top to bottom we present the CP discovery potential for the three 



We have also checked that for 6*13 > 3° our results do not change dramatically. 
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FIG. 8. The same as Fig. [3 but with 6 = 180°. 



choices of the two NSI moduh input values, (|ee^|, |eer|) = (10~^,10~^), (10~^,10~^), and 
(10~^, 10~^), respectively. The shaded regions represent the area of the parameter space in 
which CP violation can be distinguished from CP conservation at the 99% CL (3 d.o.f.). 
The yellow (light gray) regions have been obtained for ^13 = 3° using Eq. The cyan 
(dark gray) regions, on the other hand, have been obtained after searching for intrinsic 
degeneracies in ^13 and then marginalizing over ^13, using Eq. ([7]). 

We first discuss the case of ^13 = 3° (yellow regions). The results obtained are quite 
similar for 6 = and 180°. In both cases the CP discovery potential for small NSI parameters 
corresponding to {\ee^\, Icerl) = (10~^, 10~^) (top row) vanishes for the IDS25 and the 1B50 
setups, and it is non- vanishing only for two very small regions at \4>etJ.\ — 90°, independently 
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of the value of 0er, for the IDS50. The IDS25 setup presents no CP discovery potential at 
all for the case where (|ee^|, |eer|) = (10~^, 10~^) either. For this choice of NSI parameters, 
on the other hand, the IDS50 and the 1B50 setups show different performances: the former 
is able to discover CP violation due to NSI when both NSI phases are nearly maximal, 
(|0e/i|, l^erl) ~ (90°, 90°); the latter can estabhsh NSI-induced CP violation for l^erl ~ 90° 
(regardless of the value of 0e^t)- In the case of "large" NSI parameters, (|ee^|, |eer|) = 
(10~^, 10~^), the three setups yield a very good CP discovery potential for ^13 = 3°. In this 
case we observe a similar pattern of the sensitivity regions for all three setups: most of the 
space is covered, apart from two strips around (pefi = and 180° whose widths mildly vary 
for differing setups. 

As it was already mentioned in Sec. lIVBt the dependence on and eer in the z/g — )■ z/^ 
oscillation probability is quite different: the coefficient of the ee/^term is roughly one order of 
magnitude larger than the corresponding coefficient of the term in the considered range of 
energy and baselines. Therefore, the behaviour of the CP discovery potential due to NSI at 
the IDS25 and IDS50 are primarily determined by eg^if |ee^| and {eerl are comparable. This 
is clearly seen by the vertical bands shown in the top and bottom panels in Figs. [7] and [HI 
the CP discovery potential presents practically no dependence at all on eer , its behaviour 
being dominated by , since for these panels both parameters are of the same order of 
magnitude (and, therefore, the latter is suppressed by its coefficient). In the asymmetric 
case, however, (|ee^|, |eer|) = (10~^, 10^^) (middle row), these vertical bands disappear since 
now both parameters are competitive. A particularly interesting feature is seen in the right 
panels, which correspond to the performance of the 1B50 setup. In this case, due to the 
presence of the silver channel, the roles of eg^and egrare interchanged: two horizontal bands 
appear at |0er| = 90°, showing that the behaviour of the CP discovery potential is dominated 
by egT-in spite of the low statistics of the silver channel. This is a consequence of the enhanced 
role of the silver channel due to the asymmetric choice of and eer and by the size of their 
corresponding coefficients (inverted with respect to the golden channel). For |eg^| = 10~^ 
we see no significant difference for 5 = or 180° in the case ^13 = 3° (yellow regions). 

Now we discuss the CP discovery potential for the case of 6'i3 < 3°, depicted as the cyan 
regions in Figs. [7] and [H The first thing we notice is that, when marginalization over 613 and 
^13 is performed, the CP discovery potential is partially lost for all the setups under study. 
This effect is due to the presence of interference terms in the form ^i3|eeQ,| x expi{6 + (pea) 
in the probabilities. Notice that this effect is much worse in the case of 5 = 180° than for 
6 = 0. In the former case, only small islands in the {(pefi, ^'er) space survive at around the 
maximally CP-violating values of the NSI phases only for the 1B50 setup, while the IDS25 
and the IDS50 setups show no discovery potential at all for any of the considered input 
values of the NSI moduli. The marked difference between the shaded regions for 5 = and 
5 = 180° illuminates very well how complicated the interplay among the three CP-phases 
is; once we marginalize over ^13, flipping the sign of e*"^ leads to a cancellation between the 
standard and non-standard CP-violating contribution, which results in a heavy loss of the 
CP discovery potential of the facilities. This cancellation is less effective when the 1B50 
setup is considered, as the silver channel is enhanced when the golden channel gets depleted 
and viceversa. 

It is remarkable to see that high enough neutrino energies turn out to be of key importance 
in order to observe NSI-induced CP violation; the IDS50 setup always performs better than 
the IDS25, for all the input values we have considered for the NSI moduli. 

To conclude this section, we point out the following two features: firstly, as this analysis 
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has been performed with fixed mass hierarchy, it is expected that the CP discovery potential 
could become worse when the sign-Am^ degeneracies are taken into accoun10- Secondly, a 
more refined analysis of the correlations between the three phases with a proper treatment 
of backgrounds, systematic errors and marginalization over atmospheric parameters should 
be performed to confirm robustness of the observed features. 



B. 6cp fraction: Non-standard CP violation in presence of z^SM CP violation 

In the previous section, we have presented the two dimensional slice of the three- 
dimensional "CP sensitivity volume" at the very particular points 5 = or 180°. For 
different values of 6, the CP discovery potential changes dramatically due to the correla- 
tions between the three phases. We have indeed found very different features depending 
on the considered setup, the choice of the input values for the NSI moduli {\ee^\, |eer|), 
and the mixing angle ^13; in particular, many "holes" appear indicating regions inside the 
three-dimensional parameter space for which we are not able to distinguish a CP-violating 
input from the CP-conserving points. The existence and position of these holes change for 
the three setups and for the considered choices of (|ee/i|, l^eA) and ^13. Therefore, to repeat 
the procedure adopted in Figs. [7] and [8] and draw infinitely many slices of the CP discovery 
potential for different values of 6 would neither be practical nor shed any useful light over 
the intimate structure of the correlations. 

We therefore introduce a new quantity, the '^Sqp -fraction contour in the (^e^i, 0er) space", 
to condense the information. This quantity, denoted as Fs{(f)efi, 'Per)-, is defined as the fraction 
of possible values of 5 which fall into the region where CP violation can be established at the 
99% CL (3 d.o.f.) for a certain point in the (0e/i, ^er) space. Notice that this is nothing but 
the usual CP- fraction redefined on the two-dimensional plan^ (^e/x ;0er )• Such contours 
are shown for ^13 = 3° in Fig. [9l using Eq. ([6]), and in Fig. [TOl using Eq. ([7]). 

In Fig. m the white, yellow, cyan, pink and red regions correspond to Fs{(pefi, <PeT) < 
40%, 60%, 80%, 90% and 95%, respectively, for fixed ^13 = ^13 = 3°. Globally, the IDS50 
setup has the best CP discovery potential, while the IDS25 and 1B50 setups yield compara- 
ble results (the 1B50 performance being slightly better). The top and bottom panels, which 
correspond to (|ee^|, |eer|) = (10~^, 10~^) and (10~^,10~^), respectively, display similar fea- 
tures as those we have seen for 5 = or 180°. Fs{(f)e^l, (per) larger than 40% (60%) is achieved 
in almost the whole plane for all the setups in the case of "small" ( "large" ) NSI parameters. 
A Fs{(j)efj,, (per) greater than 80% is achieved for all the setups under study around |0e^| ~ 90° 
for "large" NSI input values. 

On the other hand, the result for asymmetric NSI moduh, (|ee;x|, |eer|) = (10~^, 10~^), is 
peculiar. We can see that the CP discovery potential for the IDS25 and the 1B50 setups 
is worse than for "small" NSI moduli, since a larger white region with -^^(^e/x; ^er) < 40% 
arises for ~ —90°. At the same time, the discovery potential for the IDS50 setup is 
maximal at (p^r ~ 90°. The fact that the main features of Fs{(f)efi,4'eT) are determined by 
(per rather than by appears to be again a consequence of the asymmetric choice of NSI 
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Notice, however, that in the two-baselines setups for measuring the j/SM parameters the magic basehne 
detector is able to solve the sign degeneracy in most of the parameter space It is an intriguing 

question to examine to what extent it continues to hold with NSI. This question has been investigated in 



42| but only partially and needs to be examined further 
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The concept of CP-fraction was introduced in Refs. 80|, |8l| to compare in a condensed form the per- 



formances of different proposals regarding the measurement of a given observable. It is defined as the 

fraction of the 5-parameter space (i.e., the fraction of 2n) for which a given setup is able to perform a 

94 

given task. 
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FIG. 9. Contours for the (5cp-fraction Fs in the {(j)eii, (per) plane (as defined in the text) for ^13 = 3°, 
obtained using Eq. ([6]). From left to right: results for the IDS25, IDS50 and 1B50 setups. From 
top to bottom: {\ee^^\, \eer\) = (10^^, 10"^), (10^^, 10"^) and (10^^ lO"^). The white, yellow, cyan, 
pink and red regions correspond to Fs{(pefi, (per) < 40%, 60%, 80%, 90% and 95%, respectively. 



parameters. Having \eer\ an order of magnitude larger than |ee^|, the two parameters play 
an equally important role and huge correlations arise between them, which could explain 
the absence of a peak in sensitivity at (p^fi — 0er — —90° for the IDS50 setup. When the 
two parameters are of the same order, always dominates over e^r, in agreement with our 
previous results. A comparison between left panels in Fig. [9] and, for instance. Fig. 12 of 
Ref. H indicates that marginalization over NSI hurts the CP sensitivity in a significant 
way. 
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FIG. 10. Contours for the (5cp-fraction Fg in the {(j)e^,4>eT) plane (as defined in the text) after 
searching for intrinsic degeneracies in ^13 and marginahzing over ^13, using Eq. ([7]). Left (right) 
panels: results for the IDS50 (1B50) setup. Top panels: (|e 



e/ili \ '^er | 



(10 ^,10 ^); bottom panels: 

(|ee^|, I fieri) = (10~^, 10~^). The white, light yellow, yellow, cyan, pink and red regions correspond 
to Fs{(pe,„<PeT) < 5%, 40%, 60%, 80%, 90% and 95%, respectively. 



In Fig. Uniwe show contours corresponding to Fs{(/)e^, (per) < 5%, 40%, 60%, 80% and 95% 
in white, light yellow, yellow, cyan, pink and red, respectively, after marginalization over ^13 
and ^13 below 3°, using Eq. ([7]). Results for the IDS25 setup are not presented because we 
have found vanishing Fs for any choice of \eefj,\ and |eer| under consideration. Similarly, we 
do not show any results corresponding to (|ee/i|, |eer|) = (10~^, 10"'^) either, since we have 
found vanishing Fs for all the three setups. Notice that we have included in this figure an 
additional contour for Fs < 5% in order to achieve a better resolution in small Fs regions. 
When we compare Fig. HU] to Fig. [H] we observe some contrived features. That is, the CP- 
fraction obtained for the ISD50 is not always better than that for 1B50. In certain limited 
regions the CP-fraction is actually larger for the 1B50 setup, but the area covered by the 
regions where some sensitivity to CP violation is achieved (light yellow) is larger for the 
IDS50 setup. It appears that this feature is an outcome of the complicated correlations 
between these phases. 



26 



C. CP volume fraction 



If we focus on the most general possible case in which CP violation comes both from the 
uSM as well as from NSI, it would be interesting to understand its global features. For this 
purpose we define a new quantity which we call the "CP volume fraction" . This is defined 
as the fraction of volume in which CP-violating signal can be distinguished, at a given CL, 
from a CP-conserving one in the three-dimensional parameter space spanned by 6, (pe^ and 
0eT- We use 99% CL (3 d.o.f.) to define the CP volume fraction presented in the figures in 
this subsection. 




[I, 
I 

u 




FIG. 11. Comparison of the CP volume fractions (as defined in the text) as a function of ^13, 
using Eq. ([6]). From left to right: IDS25, IDS50 and 1B50 setups. Results for (|ee;,|, |eer|) = 
(10~^, 10"^), (10~^, 10~^) and (10"^, 10^^) are represented by dotted, dashed and solid lines, re- 
spectively. 



In Fig. [m we present the results of the CP volume fraction computed using Eq. 
for the three setups and the three choices ("small", "asymmetric" and "large", as usual) 
as a function of ^13. Consider first the IDS25 and IDS50 setups: we can see that, at both 
setups, the CP volume fraction achievable for (|ee^|, Icerl) = (10"'^, 10~^) and (10~^, 10~^) are 
extremely similar. This means that the CP discovery potential for both setups is dominated 
by |ee^| and the effect of |eer| is marginal, in agreement with the results obtained in previous 
sections. On the other hand, we see that when \eefj,\ = 10~^ (solid line), the CP volume 
fraction rises abruptly to approximately 80% or above for ^13 > 2°. 

This is not the case for the 1B50 setup: we can see that in the case of "small" NSI 
parameters (dotted line) it shows a similar (but slightly worse) sensitivity than that of the 
IDS50. The CP volume fraction is approximately 50% for ^13 > 2°. However, as soon as 
leerl = 10~^, it rises to 70%, and reaches 90% for larger |ee^| = 10~^ at larger 6^13. The 
behaviour clearly indicates that the 1B50 setup is taking advantage of the combination of 
both channels in an efficient way. 

Another interesting observation that we can draw from Fig. [11] concerns the 613- 
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dependence of the CP volume fraction. Firstly, the CP volume for ^13 = reflects the 
ability of each setup of observing CP violation exclusively due to NSI. While for the IDS25 
setup the CP volume is nonzero only when |ee^| is "large", the higher energy setups would 
be able to measure a CP violating signal even in the case of "small" NSI parameters. This 
confirms the results presented in Sec. IVI Al (Figs. [7] and |H])- Secondly, we can see in the 
"small" NSI parameters case, and also in the "asymmetric" case for the IDS25 and IDS50 
setups, that the CP volume fraction decreases at "large" ^13, reaching a maximum for 
613 G [1°,2°]. It means that, when the NSI parameters are rather small, correlations with 
a "large" ^13 can actually reduce the CP-discovery potential of those facilities. This is no 
longer true when the NSI parameters are set to be "large": in this case, the CP volume 
fraction increases for increasing 6is value but remains practically unchanged for ^13 > 2°. 
No destructive correlations arise when the NSI parameters are sufficiently "large". 




FIG. 12. Comparison of the CP volume fractions (as defined in the text) as a function of ^13 
for {\€e^,\,\eer\) = (10-3,10-2). From left to right: IDS25, IDS50 and 1B50 setups. Results 
for the first detector (MIND at the Intermediate Baseline) and for the second detector (either 
MIND at the Magic Baseline or MECC at the Intermediate Baseline) are shown in red and green, 
respectively. Combination of the two detectors (Vi+2) see text) is shown in blue. The dotted 
black line corresponds to the simple sum of the CP volume fractions of the separate baselines or 
detectors, Vi^2 (see text). The dashed black lines represent the synergy, i.e. V1+2 — 14,2 (see text). 

To conclude this section, we address here the issue of the relative importance of the differ- 
ent detectors and different baselines, and of the possible synergies between them. We present 
in Fig. [T2]the CP volume fractions computed for each detector and their combinations for the 
three setups. Only the case of "asymmetric" NSI parameters, (|ee^|, |eer|) = (10"^, 10"^), is 
shown. Notice that it is precisely the case where the synergy between channels/detectors is 
most important because the two NSI parameters are equally relevant in the golden channel, 
which presents the largest statistics. 

Let Vi be the CP volume fraction obtained by analyzing data given by a detector i, where 
i = 1 refers to the MIND detector at L = 4000 km and z = 2 to the MIND detector at 
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L = 7500 km (left and middle panels) or to the MECC detector at L = 4000 km (right 
panel), respectively. The volume Vi is represented by red lines, V2 by green lines in all 
panels. The CP volume fraction obtained combining two detectors, V1+2, is represented 
by the blue lines. Notice that V1+2 is not the simple sum of the CP volume fractions of 
the two detectors, Vi + V2. When we combine the CP volume fractions of two detectors, 
we must take into account that some part of the parameter space can be covered by both 
detectors at the same time. We call the latter volume fraction as Kveriap, the fraction of 
the three-dimensional parameter space that is covered simultaneously by both detectors. 
Then the correct definition of the simple sum of the CP volume fractions of two detectors is 
Vi,2 = Vi + V2 — Overlap- The volume Vi,2 is shown by the black dotted lines in Fig. [T2j We 
can see that, in general, V1.2 does not coincide with V1+2, because the combination of the 
data from two detectors could in principle cover regions of the parameter space that are not 
accessible to each detector separately. This is nothing but the effect of the synergy between 
two detectors, an increase in the CP volume fraction due to simultaneous analysis of the 
data sets of two detectors, a procedure different from the simple sum of the sensitivities 
achieved by each one separately. The synergy V1+2 — Vi,2 is represented by the black dashed 
lines in Fig. [T2l 

The first thing to be noticed is that the synergy is never much larger than approximately 
10% for any of the considered setups. In the IDS25 and IDS50 setups, Vi,2 (black solid 
lines) almost coincides with Vi (red lines for the intermediate detector). This means that 
most of the CP volume fraction due to the magic baseline detector is already covered by 
the intermediate detector. In the case of 1B50 setup the situation is different. For ^13 < 1° 
the MECC detector does not contribute at all, as the red, black and blue lines coincide. 
For > 1°, the MECC contribution grows linearly and it starts to cooperate with the 
MIND detector, as expected. However, we see that Vi^2 (black) and (blue) are closer 
compared to the cases of IDS25/50. This means that the overlap region V^veriap is relatively 
small in this case. The two detectors are complementary, as they test different regions of 
the parameter space. 

VII. SUMMARY AND CONCLUSIONS 

In this paper, we have performed a complete study of the effects of NSI in neutrino 
propagation at the High Energy Neutrino Factory (HENF). We paid particular attention 
to the correlations among the whole set of oscillation parameters, the uSM and the NSI 
ones. Our analysis is the first one which takes into account all the NSI parameters at the 
same time in the simulations. Among all the new facilities proposed and discussed in the 
literature, we have chosen HENF due to its high energies and very long baselines, as it 
is likely to be best suited to study effects of NSI in propagation. We have examined the 
three setups defined in Sec. IIIIDI (a) IDS25, the standard IDS setup; (b) IDS50, a 50 GeV 
upgraded version of IDS25; and (c) 1B50, a one baseline 50 GeV setup with a composite 
detector capable of detecting z/^ as well as z/^. The comparison among the three different 
setups has been performed keeping in mind that a possible optimization of HENF could be 
performed in order to search for New Physics (NP) in case that ^13 turns out to be at reach 
by the ongoing (or forthcoming) neutrino experiments. 

Our analyses have been divided into two parts: in the first part (Sees. [IV] and |V]) the 
sensitivities to NSI parameters have been studied after a brief discussion of effects of NSI on 
the measurement of i^SM parameters; the second part (Sec. rVI|l is devoted to the analysis 
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of the CP violation discovery potential of the setups under study associated to the three 
phases, 6, (pe^ , and (per ■ 

Results for the first part can be summarized as follows: 

• Significant correlations between 6^13, and e^r have been found at HENF setups. 
Such effects can be reduced placing a detector at the magic baseline, but they cannot 
be eliminated. The ultimate sensitivity to 6'i3 is worsened by a factor 3, 5 and 10 
for the IDS50, IDS25 and 1B50 setups, respectively, when NSI are included in the 
marginalization procedure compared to the case without NSI (see Fig. [2]). On the 
other hand, no correlations with e^r and only a marginal effect of e^a have been 
observed. 

• The sensitivities to \ee^\ and |eer| essentially come from — > (and, to a lesser 
extent, z/g — )■ v^) oscillations channels, while the sensitivities to \t^r\ and the diagonal 
NSI parameters are achieved mostly through the z/^ — )■ i/^ disappearance channel. We 
have numerically checked that these two sectors are practically decoupled, since no 
significant correlations on their sensitivities have been found. The marginalization 
over does not affect significantly the sensitivity to ee^, however the converse is not 
true. 

• A higher sensitivity to \€e^\ is achieved by the higher energy setups, IDS50 and 1B50. 
For these setups, the sensitivity to |ee^| at 95% CL (2 d.o.f.'s) for ^13 = (3°) is in the 
range [0.5, 1.7] x lO^'^ ([0-7, 2.5] x 10^'^), depending upon 0e^. The sensitivity obtained 
by the IDS25 setting is worse by a factor of 3 (see Fig. [3]). 

• The sensitivity to |eer|, which is worse than the sensitivity to |ee^|, varies significantly 
between the different setups, although generally better for the two baseline setups. 
The smallest value of |eer| that can be excluded at 95% CL (2 d.o.f.'s) for ^13 = 
(3°), depending on (per, is in the range [4, 5.5] x lO'^ ([5.5, 12] x lO'^) for the IDS25, 
[2,3.3] X 10-3 ([4, 6] X 10"=^) for the IDS50 and [5,10] x 10-3([6,17] x lO'^) for the 
1B50 (see Fig.S]). 

• The sensitivities to the diagonal parameters and to e^r, on the other hand, are quite 
independent from the setup under consideration. The sensitivity to (eee — Ctt) is quite 
limited, ~ 0.2 (~ 0.1) for ^13 = 0(3°), since this parameter shows up only at 0{e'^) 
in all the oscillation channels. The sensitivity to this combination of parameters is 
strongly limited by the uncertainty on the matter density in the earth. 

The sensitivity to (e^^ — e,-^) is approximately one order of magnitude better and 
reaches ~ 0.03 (~ 0.05) for ^23 = 45° (^23 = 43° or 47°), see Fig. d The sensitivity 
to |e^r| shows a remarkable dependence on the corresponding phase (p^r (much larger 
than in the case of \ee^\ and leerD- The values that can be tested range from (9(10~^) 
to C(10-3) for \(p^r\ = 90° and \i)f,r\ = 0, 180°, respectively (see Fig. [5]). 

In the second part of this work, we have studied the potential of the HENF to discover 
CP violation due to the uSM phase and/or the NSI phases <pefj.,(peT- This has been done 
by computing the discovery potential in the three-dimensional space, i.e. the region of 
sensitivity to CP violation at the 99% CL (3 d.o.f.'s) in the {S, (pe/j,, (per) space. Since (p/j^ris 
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not expected to be correlated to the above mentioned CP-phases we have not studied the 
correlations between (p^r and them in this work0 

We have distinguished two different cases depending on the value of 6*13: either 6^13 is 
at reach in the currently running (or in the next generation) neutrino experiments, or it 
is out of reach. In the former case ^13 has been fixed to be 3° as a typical value, and in 
the latter, we have marginalized over ^13 in the allowed range, ^13 < 3°. We have chosen 
"reasonable" values for the NSI moduli, (|ee;_j|, |eer|) G [10~^, 10~^], paying special attention 
to the particular case where |eer| is an order of magnitude larger than \eefj,\- In this case 
effects of both parameters are competitive in the golden channel and interesting correlations 
arise between them. Results obtained in this section can be summarized as follows: 

• li 6 turns out to be at a CP-conserving value {i.e., when CP violation is exclusively 
due to NSI) a strong dependence on the true values of 613 and S (either or 180°) 
is observed. For ^13 = 3°, CP violation can be established in most of the parameter 
space only for "large" NSI parameters, {\ee^\, |eer|) = (10~^, 10~^) for both 6 = and 
180° in all the three setups. When marginalizing over ^13, however, CP violation can 
be established in a significant amount of the parameter space only for 6 = and 
only for the IDS50 and 1B50 setups. For 6 = 180° destructive interference occurs 
and CP violation can be distinguished from CP conservation only in the points where 

(l^e^l, I0er|) = (90°, 90°). 

• For generic values of the z/SM CP phase, the 6cp fraction Fs (defined in Sec. IVI Bp 
is larger than 60% for ^13 = 3° and "large" NSI parameters in all (most) of the 
phase space for the IDS50 and 1B50 (IDS25) setups. Yet, the regions in which 
Fg > 40% extend to most of the phase space even for "small" NSI parameters, 
(|ee/i|, Icerl) = (10"'^, 10"'^) , for all the setups. When marginalizing over ^13 the CP 
sensitivity coverage becomes much poorer, though. Vanishing CP sensitivity {Fg = 0) 
is found for the IDS25 setup for all the input values of the moduli. The same result 
is obtained for the IDS50 and the 1B50 setups for "small" NSI parameters, and non- 
vanishing Fs is found only in some regions for the "asymmetric" and "large" choices 
of the NSI moduh. 

• We have also analyzed the CP volume fraction, defined in Sec. IVI CI The relevant 
features are as follows: (1) The high energy setups have higher CP volume fraction 
than the IDS25 for all cases considered; (2) The CP volume fraction at the IDS25 and 
the IDS50 setups is determined by the value of |ee^|, whereas \eer\ has no impact; (3) 
For "small" values of the NSI parameters, the CP volume fraction shows a maximum 
at ^13 ~ 1° at all setups (the same is true also for "asymmetric" NSI in the IDS25 and 
IDS50 setups). For |eer| = 10~^ the 1B50 setup has the largest CP volume fraction 
among the three setups, ~ 70% (90%) for ^13 > 2° for \eef,\ = 10-^(10-2). 

• Finally, we have also studied the synergy between the two detectors located at different 
baselines (IDS25 and IDS50) and between the golden and silver channels (1B50) for 
"asymmetric" NSI parameters, since the synergy may be most prominent in this case. 
We found the synergy to be most significant for the IDS50 setup, though it is still 
small (roughly speaking, 10%). 

Notice that our results for the CP discovery potential constitutes only a first step toward more careful 
studies of the complicated correlations arising between S, (jie^i and (per- Analyses which include a proper 
treatment of backgrounds, systematic errors and marginalization over other parameters (with particular 
emphasis, in this case, on the sign of the mass hierarchy) should be performed to confirm the results 
presented here. 01 



To summarize the outcome of our comparison between the three setups we make the 
following remarks: (1) generally speaking, the high energy setups, IDS50 and 1B50, are 
better than IDS25, a naturally expected result since NSI behave as a generalized matter 
effect; (2) the former two settings have their own merit and demerit: their sensitivities to |ee^| 
are comparable, while the one to |eer| is higher for the IDS50. An accurate measurement of 
^13 is more robust to a possible obstruction by NSI for the IDS50. The CP violation discovery 
potential is higher for the 1B50 setup, in particular for large ^13 ~ 5°. However, the IDS50 
setup may yield a better global performance than the 1B50 because of its robustness thanks 
to the detector located at the magic baseline. 

As a final remark, it should be noted that the sensitivities achieved at the high energy 
setups studied in this paper are remarkable, and close to the edge of the effects produced by 
some neutrino models of New Physics. If New Physics is at the TeV scale, it is quite likely 
that new sources of CP violation beyond that of z/SM exist. If one interprets our exercise as 
an example for such generic cases, apparently, the lessons we have learned are that interplay 
between these phases are highly nontrivial, and any regularities between them or knowledge 
of the right model would be of crucial importance. What we have definitely learnt in this 
study is that higher neutrino energies (such as 50 GeV) have proven to be crucial in order 
to pursue these elusive NSI effects in neutrino oscillation experiments. 



NOTE ADDED 

After submission of this paper, an excess of z^appearance events far above background 



was reported by the T2K [82| and the MINOS [83(] experiments, indicating a large value of 
613. A global fit which includes these new data gives ^^13 ^ 8.3° 8^. Therefore, we have 
performed a first scrutiny of our results under the premises of such a large value of 6*13, which 
is larger than what we examined during our analysis, i.e., ^13 up to 5°. The main outcomes 
of this preliminary analysis are the following: 



1. the maximal sensitivity to ee/j. and eer is not affected significantly at any of the con- 
sidered setups although the results are, in general, slightly worsened; in addition, the 
sensitivity minima are slightly displaced in the (per axis with respect to the 613 = 3° 
case; 



2. the sensitivity to e^^ and to the diagonal parameters is not affected significantly, 
either; the only exception being the sensitivity to (e^^j — e^r) at the 1B50 setup, which 
gets worse by a factor ~ 3; 



3. from a very preliminary scan of the CP discovery potential, we get slightly worse 
results since a large 6^13 gives stronger correlations between the NSI phases and 6, (as 
mentioned above). These preliminary results are in reasonable agreement with the 
extrapolation of Fig. [T21 



Since full understanding of effects of NSI in propagation in the case of large ^13 requires 
new theoretical machinery and the analysis can go beyond what we have presented in this 
paper, we leave the detailed discussion of this case to a future communication. 
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Appendix A: Expressions of Expanded oscillation probabilities 

In this Appendix, we show the oscillation probabilities Pq,^ in matter with constant den- 
sity, in presence of NSI affecting onlyto propagation in matter. We start from the oscillation 
probability expansions derived in |67| where e^/j, 6*13, and [S.m\^l lS.m^^ are considered the 
expansion parameter^. Here we expand also on bQ2?, = ^23 — 7r/4, considering therefore: 

^ ■ {ea/3, di3, Am^i/Am^^, 5023} 

as the order e expansion parameters. 

The oscillation probabilities for golden and silver channels at O(e^) are given by the 



Remember that, in the set-up considered in this paper, A ^ A31. 
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following formulae: 
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where A^^ stands for the standard oscillation amplitude, A for the matter density, Ajj 
(Am? /2£'), and Sij and Cij stand for sin% and cos%, respectively. 
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On the other hand, since the sensitivities to e^^ and e^T are mainly achieved through the 
golden and silver channels, for the z/^ — sector we show here only the dependence on eaa. 
and e^T- of the relevant oscillation probabilities, which will be called P^n^'- 



yNSI _ pNSI 



tJ.T 



{^3 ( ) + Re{e^,)}{AL) sin (AgiL) 



A 
2 A 



sm 



- (Re(e^.))^ (Aiy cos (A31L) - (Im(e^.))^ l^i^L) sin (A31L) . + 0(6 



A 
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(A3) 



The complete oscillation probabilities at quadratic order in e can be found in 6^ ■ 
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